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In troduction 

To assist investigators in locating the “critical points" (points 
of highest magnitudes) in rotating components whole field techniques 
for strain (or stress) analyses are needed. One such technique is a la- 
ser speckle technique that produces Young's fringes. The extension of 
this speckle technique to the strain (stress) analysis of rotating parts 
is the object of this paper. 

In extending this speckle technique to the analysis of rotating 
parts use is made of a film plate "sandwich" procedure [1,2]* from which 
certain advantages are realized. One advantage over point techniques 
is that the sandwich procedure retains the full field information found 
in usual speck! egrams . An advantage relative to holography is that the 
information available from specklegrams allows direct evaluation of in- 
plane displacement which is needed to calculate in-plane stresses [3,4]. 
An advantage of sandwich specklegrams over usual specklegrams is that 
the range of measurement can be increased to include displacements on 
the order of a few microinches up to values approaching 0.1 inch. Also, 
positive or negative displacement can be assigned to a point using the 
sandwich speckle procedure, something that usual (single film) speckle- 
grams cannot reveal. Timing constraints for gathering displacement data 
from a rotating component are considerably relaxed with the sandwich 
speckle interferometric procedure as compared to holographic and conven- 
tional speckle techniques [5,6]. 


* Numbers in brackets refer to references . 
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Finally the analysis of the data for determining strains is more 
convenient than holographic analyses. 

Essentially the speckle procedure involves photographing the speck! 
produced by a pulsed laser when the laser light impinges on a rotating 
structural component (or part). In this sandwich procedure the compo- 
nent is then stopped and repositioned as close as one can in the same lo 
cation that it occupied at the time of the first photograph. A second 
photograph is then made using the pulsed laser to illuminate the station 
ary part. The two negatives containing the speckle patterns are sand- 
wiched together with the emulsions in contact and the images realigned 
relative to each other so that Young's fringes can be observed when a 
small area of the sandwich is illuminated with a continuous laser. Dis- 
placement components can be calculated from the Young's fringe data. 
Stresses can be calculated by differentiating the displacement data. 

Analytical Description of the Sandwich Speckle Fringes : 

Consider a point, 0, contained in an aperture and illuminated by a 

coherent light source, S, as illustrated in Fig. 1. The light wave 

scattered by 0 to a point P on a distant screen maybe represented 

by the real part of 

__ _ i(wt +$) 

dE po = Ae dxdy (1) 

A = amplitude of the optical disturbance, 
w = transverse frequency of light, and 
$ = phase angle of the optical disturbance. 

Let point Q be a nearby point located relative to 0 by vector 8 
B = xi + yj + zk 


where 


( 2 ) 
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The optical disturbances 0 and Q are out of phase by at 
point P. The phase difference, Acj> , is proportional to the difference 
in path lengths along SOP and SQP and is represented as 

A<J> = — [(QS • QS) 1/2 + (QP • Qp] /2 - (OS • 0SI /2 - (OP • Op}/ 2 ] (3) 

where 


x = the wave length of the coherent light source. 

It is noted that 

QS = 03 - 3 and (4) 

QP = OP - B (5) 

By substituting Eqs. (4) and (5) into Eq. (3), using a bionomial 
expansion and neglecting higher order terms (assuming |ET|«| OP | , (OS’!) 


n 

we obtain A^-r- 2 - 


OS 


OP 


(OS • 0S) 1/2 (OP • 0P) 1/2 J 


B 


n s = a unit vector along OS = - k 


n Q = a unit vector along OP = ii + mj + nk 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

where i, m, n are direction cosines of n Q . The optical disturbance at 
P due to Q is now expressed as 


taking 

and 

then 

or 


2tt 


a* =- - (n s + n Q ). B 


At » - ~ [xz + ym + z(n-l)] 


dE PQ 3 Ae ’ 


i jwt + t ~ [xa + 


ym + z(n-l)] 


j-dxdy 


(ID 


The combined optical disturbance at P from all points in a two dimen- 
sional aperture is 
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p AG£ is 

quality 


E,Ae 1 '< wt+ * ) fje’ 1 T [ “ + my + z(;n - 1);l dxdy (12) 

Consider a small aperture ('say aperture 1) centered at the origin and 
contained in the X - Y Plane. The disturbance at P due to aperture 1 
is 

E. - 7te i(wt+ + ) / e" 1 + * m] dxdy (13) 

area 

Next consider a second aperture the same size as aperture 1 displaced 
relative to the first by ax. Ay, az. The disturbance at P due to this 
aperture is 

£ = f e ‘i ^ Cfx+ax) JL + (y+ay) m + Cn-l)az] dxdy 

~ 2 ' area 


or since ax. Ay, az are constants 

E 2 = E 1 e” 1 * 'T ^ xz + Ayrn + A2 ^ n “ 1 ^ 

■T-he total disturbance at P then is 

T - 7 ] l + p* 1 H + fiyn> + <iz(n-l)] ) 
t T0TAL “ t l l 1 e * / 

The intensity at P is given by 

! P = C E TOTAL ' E *T0TAL 


(15) 

(16) 


( 17 ) 


where 

C = proportionality constant 

1e __ 

E T0TAL = complex conjugate of E TQTAL 
Expanding Eq. 17 and collecting terms yields 


(1C) 

where I is the intensity due to a single aperture and typically repre- 
1 

sents a halo while the terms in parentheses represents "Youngs fringes" 


Ip = 21^ / 1 + cos [ZAX tmAy + (n - 1) Az] 
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across the halo. 

For further consideration of the sandwich technique let the X-axis 
be coincident with the displacement component being measured (jnormal to 
the camera axis) and let Az represent the "gap" between film plates along 
the line of sight of the camera. 

Then 

I p = 2^ 1 1 + cos [mx + (n - 1)az] \ < 19 > 

Fringes occur when 

MX + (n-l)AZ - /- 2Q s 


with N = + 1 , + 3, - - - 

If az - 0 or if ax »a z then Eq. 20 reduces to 

A z% (2i) 


or 


AX 


Nx(x 2 + z 2 ) 1/2 
~ 2 x 


( 22 ) 


and is equivalent to the usual double exposed specklegram on a single' 
film plate (Fig 2). 

If ax =0, Eq. 20 reduces to 


(n-1 )az 


N_X 

2 


(23) 


which describes a fringe pattern of concentric circles (Fig 3). The 
spacing between these circles is inversely proportional to az. 

If ax and az are comparable in magnitude, the resulting fringe pattern 
is curved (Fig 4). 

Data can be interpreted according to Eqs. (20), (22) or (23). If 
the fringes are relatively straight Eq. (22) yields the desired informa- 
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Figure 2. Interference Pattern From Displaced Speckle Patterns 




8 


ORIGINAL PAGE tS 
OF PO >R QUALITY 



Figure 3. Interference Pattern From Separated Speckle Patterns 
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Figure 4. Interference Patterr. From Displaced and 
Separated Speckle Patterns 
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tion, If the curvature of the fringes is judged to be significant, 
translate and/or rotate one speckle photograph with respect to the other 
until the fringes become concentric circles and evaluate az from the 
concentric rings and Eq. (23). It is noted that n = z/(x 2 + z 2 }^ 2 . 

Now with the determined value of az and Eq. (20) ax can be evaluated. 
Alternatively, if the curvature is significant, translate the plates to 
produce closely spaced Youngs fringes then measure the fringe spacing 
between adjacent fringes at the center of the halo where (n-1) stays 
approximately zero. Hence az need not be determined explicitly. By 
letting the second specklegram (the one furthest away from the light 
source) correspond to a particular specklegram, say the one of the load- 
ed model, more information is available. In this case az is positive. 

Note that V, (which is X/R) is positive on the positive side of the 
X axis and negative on the negative side of the X axis, which means 
that fringes are positive on the positive side of the X-axis and negative 
on the negative side of the X-axis (see Fig. 5a). Fig 5(c) represents a 
typical curved fringe pattern. Figs. 5(a) and 5(b) illustrate that the 
resulting pattern Fig. 5(c) is composed of two parts. One part repre- 
sented by Fig. 5(a) is due to a displacement along the x axis. The other 
part (Fig. 5(b)) is due to a seperation of corresponding speckles (or a 
seperation of the two emulsion film plates.) When part (a) is added 
to part (b) a curved fringe pattern results as shown in part (c) and 
indicates that the displacement is toward the center of curvature of 
the fringes. Part (d) illustrates a negative displacement and when add- 
ed to the concentric circles of (e) caused by a positive az results in 
curved fringes illustrated in (f). Again note that the displacement 
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\s towards the center of curvature. Therefore, it is evident that curved 
fringes can be used to indicate the sign of the displacements. If a 
fringe pattern is not curved when using the sandwich technique and the 
positive or negative nature of the displacement at a point is needed az 
can be induced to reveal this information. 

Next there are cases when displacements are too small to cause fring- 
es in a halo. For example in a cantilever beam the displacement near the 
fixed end is small and no Youngs fringes will be evident in the halo cir- 
cle. In this case the region near the fixed end is the region containing 
the highest strains and is thus the most important region from a design 
view point. This difficulty can be alieviated with application of the 
sandwich technique by mechanically displacing one film plat* relative to 
the other to produce more closely spaced Youngs fringes. The additional 
fringes correspond to an apparent rigid body movement and do not contri- 
bute to the calculated strain components. Normal strain components are 
expected to be more accurately evaluated than shear strains since rigid 
body rotations do alter the values of the cross derivatives of displace- 
ments. Therefore to evaluate the strain state at a point ? it is recommend- 
ed that the normal strain be determined in three directions, say 0°, 45°, 
90°, and these values used to calculate shear strain as is done with strain 
gage rosette calculations. 

Experimental Data Collection : 

To explain how data for rotating structural members can be collected 
via the sandwich speckle technique an example problem is presented. 

The example problem selected was made up of two rotating tensile bars 
with weights on the free end of each specimen for load and balance, Fig. 

6 . However, of more importance theoretical strain results could be 
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compared with experimental results for evaluation of this experimental 
technique. Some effort was expended to isolate the system from vibra- 
tions by using rubber pads and mounting the rotor shaft in air bearings. 

It should be emphasized, however, that the speckle technique does not 
require the same stringent vibration isolation as does holographic 
techniques. The bar assembly rotated at about 3450 rpm. 

The arrangement used to time the events for obtaining data is shown 
in Fig. 7. The most important component in thi's arrangement was the 
pulse ruby laser. First it produced coherent light necessary for a speck- 
le photograph. Secondly by being pulsed a dynamic event could be re- 
corded on film without blur. Data was gathered by taking a speckle 
photograph of one of the tensile bars while it was rotating. Next the 
tensile bar was stopped and repositioned in the same location that it 
had when the dynamic photograph was recorded and a speckle photograph 
taken. The two negatives were sandwiched together and adjusted relative 
to each other until Young's fringes could be seen as a laser light was 
passed through a small area of the negative. When Young's fringes could 
be seen the negatives were judged to be correlated. 

The sequence of events for making the speck! egrams is noted. The 
motor for rotating the tensile bars was turned on and ran for several 
minutes to obtain a stable rotation. The timing sequence for the dynam- 
ic speck! egram was initiated by pressing the single beam reset lever of 
oscilloscope I. One of the tensile bars had a small mirror glued to it. 
When the mirror reflected the CW argon laser beam (after reset) on to 
photodetector I an electrical pulse initiated a timing signal for A- 
gate (a 20 volt pulse) of oscilloscope I. The timing signal was set 
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Schematic of Events 
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so that the desired tensile bar would rotate into proper position as 
the ruby laser pulsed to make a specklegram. Time (between 1.0 and 
1.5 milliseconds) had to be allowed for the ruby laser power supply to 
energize the ruby laser. During this interval between 1.0 and 1.5 milli- 
seconds the pulse ruby laser was Q-switched and the ruby laser delivered 
a large burst of light for about 50 nonoseconds duration. The sequence 
occured as follows. The timed 20 volt signal (A-gate) from oscilloscope 
I triggered oscilloscope II. At this instant B-gate (20 volt pulse) of 
oscilloscope II turned on the ruby laser power supply. One millisecond 
after B-gate activated the ruby power supply, A-gate of oscilloscope II 
triggered the pockel cell pulse generator. The pockel cell pulse genera- 
tor sent a 100 volt switching signal (1.5 microseconds in duration) **o 
the pockel cell power supply. The pockel cell was switched off and then 
on, allowing the ruby laser to pulse. The coherent pulsed laser illumina- 
ted the rotating tensile bar and a dynamic specklegram was recorded. 

The dynamic specklegram film plate was positioned with the emulsion 
toward the camera lens and directly against a blank glass plate of equal 
thickness. The blank glass plate was used so that the light path of the 
dynamic specklegram would approximate the light path of the second (or 
stationary) specklegram, where the emulsion is faced away from the camera 
with the glass film backing in between. These steps made it possible for 
the two emulsions from the speck! egrams to be placed in direct contact 
during the sandwich process. After developing and drying the dynamic 
specklegram, it and the glass h.lank were repositioned in the film holder. 

The film holder is shown in Fig. 8. The screws located the surface of the 
glass blank that was closest to the camera lens. The glass, in turn located 


17 


ORIGINAL PAGE IS 
OF POOR QUALITY 



rew 


Screw 


Screw 


Pin 


Figure 8. Film Holder 
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the film emulsion plane. The pins positioned the glass blank and the 
film edges. The rotating bar which was photographed in the dynamic 
specklegram was repositioned in the same imaged location by viewing 
through the repositioned film. With the bar in this imaged position the 
glass blank and dynamic specklegram were removed from the holder. Now 
a stationary specklegram was taken of the stationary bar by liqht scat- 
tering from the bar, through the lens, through the film glass backing and 
activating the emulsion. This procedure as noted before made it pos- 
sible to sandwich the two speck! egrams (the dynamic and stationary speckle- 
grams) together with the emulsions touching. It should be noted that there 
was no opaque backing on the film glass of the stationary specklegram 
so that the pulsed light could pass through the glass backing and acti- 
vate the emulsion. 

Next in the sequence after developing and drying the stationary 
specklegram, the dynamic and stationary specklegrams were sandwiched 
together and held with the positioner shown in Fig. 9. This positioner 
holds one film fixed while the second film is translated and/or rotated 
relative to the first such that Young's fringes could be observed. Dur- 
ing this trial and error process if Young's fringes are observed the 
negatives are optically correlated such that the wave front generated 
through the first negative interferes with the wave front generated 
through the second negative. A typical fringe pattern is seen in Fig. 

10. When fringes were observed the two films were locked together and 
data was taken at the desired points using the scheme shown in Fig. 11. 

In this illustration with the film emulsions together the spacing, 
aZ, between the emulsions occured only due to slight variation in flat- ' 


19 


ORJQfNAL PAGE IS 
POOR QUALITY 



Figure 9. Film Positioner 
(A) 
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ness of the film plates. The amount of variation, &Z, was not needed in 
the calculations and was not determined. Sometimes the dynamic and sta- 
tionary specklegrams would not correlate. In these instances the speci- 
men (bar) 'position was rotated slightly (say 4. or 5 thousandths of an 
inch at the bar tip) and a second stationary specklegram recorded. If 
the dynamic specklegram and the second stationary specklegram would 
still not correlate a third stationary specklegram was taken in the other 
direction relative to the original position. Two to three trials usu- 
ally produced correlation. It was noted that in a controlled experi- 
ment with a continuous laser as much as 0-07 inch shift between speckle- 
grams was possible without loss of correlation. 

Strains were calculated at various points on the centrifugally 
loaded bar and compared with theoretical values. In this example, calcu- 
lations were simplified since displacements at the center of the bar 
were axial. Thus it was only necessary to find the displacement values 
at these points and to numerically differentiate to obtain axial strains. 

A second order central difference formula was used to evaluate the deriva- 
tives. We chose the data points to be 0.0625 inch to each side of the 

point of interest. The theoretical calculations were made assuming no 
* 

significant bending and treating each cross section to be in uniform 
tension. The tension was calculated using the inertia vector caused by 
the material mass from the point of interest radially outward. Results 
of three example points are tabulated in Table I. 
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Table I. Experimental Results 


Point 

Distance from 
centerline to 
interest point, 
inches 

Calculated 
• strains from 
data 
in/ in 

Calculated 
strains from 
theory 
in/ in 

%Difference Based 
on Experiment 

1 

2.125 

305 x 10" 6 

295 x 10' 6 

3.3 

2 

2.417 

289 x 10“ 6 

293 x IQ' 6 

1.4 

3 

2.750 

274 x 10' 6 

288 x 10' 6 

5.1 


The results seem 

to compliment 

each other and 

suggest the credi- 


bill ty of this approach for analyzing strains (or stresses) in rotat- 
ing parts. 

Stresses which are induced in structures because of the rotation 
of the structure can be evaluated with speckle techniques. The data in 
comparison with holography data is easy to obtain and especially easy 
to analyze. Motion or vibration of the structure is not as restricted 
with speckle techniques as with holographic techniques. The data collect- 
ed is essentially full field. 

The following paragraphs deal with numerical techniques for inter- 


preting laser speckle data. 
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Integration of Laser Speckle and 

Finite Element Techniques of Stress Analysis : 

Contemporary methods of stress analysis for design of complex struc- 
tural components generally fall into two categories: (1) Experimental 
methods, (2) Numerical techniques. Examples of current experimental 
methods include laser related techniques (holography, speckle interferom- 
etry), grid methods (Moire analysis) and photoelastici'ty . Some of the 
currently used numerical techniques include finite element modeling, 
finite difference analysis and numerical boundary integral methods. 

Both types of approaches exhibit distinct strengths and weaknesses. For 
example, the experimental methods of holography and speckle interferome- 
try yield measurements of surface displacements of a structure. The 
determination of stress typically requires differentiation of this data 
and can result in a substantial loss of accuracy. Also regions of maxi- 
mum stress are sometimes regions of relatively little displacement; thus, 
experimental data may be sparce in those areas. Numerical analysis (eg. 
finite element modeling) however, is often capable of accurate stress 
evaluation but requires an accurate knowledge of boundary conditions 
(boundary displacements and/or tractions). Such boundary conditions 
are not always known (eg. surface pressure loads on a rotating turbine 
blade). The potential exists for combining these techniques in a com- 
plimentary fashion to make more effective use of their respective strengths 
and to overcome their separate limitations. By developing procedures 
for effectively intermeshing the two types of analysis techniques a po- 
tentially more effective tool for stress analysis and structural compo- 
nent design can be expected. Investigation of this possibility consti- 
tutes the problem examined herein. 
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Related Research For Use of Laser Speckle 
and Numerical Techniques : 

Relatively little research has been devoted to developing an inte- 
grated approach to stress analysis of the type proposed herein. Tradi- 
tionally, the experimental and numerical techniques have been used as 
separate tools with application to a common problem primarily for compar- 
ison and validation of results. Occasional reports of the use of exper- 
imental data for input to a numerical model have appeared (7,8) but 
without particular effort toward developing or exploring the potential 
of this concept. Rowlands (9) reported the use of a finite element 
type discretization and curve fitting technique on experimentally ob- 
tained (Moire) displacement data for strain analysis. KIs results were 
encouraging but his method did not fully exploit the strengths of the 
finite element method in that all model displacements were specified 
from experimental data. Thus, equilibrium and structure geometry played 
no role in the strain calculations. Additionally, Rowland's approach re- 
quires a large amount of experimental data in the regions of interest. 
This requirement cannot always be met when using displacement measure- 
ment techniques. 

Recent reports of the combined use of experimental laser speckle 
and numerical boundary integral techniques have appeared, (10,11) In 
this approach a relatively small region of interest on the structure is 
"Isolated" for analysis. Experimental displacement data is collected 
around the boundary of this region and used as input data to a numerical 
model of this sub-region. Good results have been reported for the test 
problems that were investigated, A limitation of this approach is the 
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inherent constraints on the boundary integral method as a numerical 
technique. It is best suited to modeling small regions over which grad- 
ients do not change appreciably. An entire complex structure cannot be 
effectively modeled with this technique. Thus regions of critical stress 
must be known apriori. Again, a large amount of experimental data must 
be available in the region of interest. Finally, the technique is not 
easily adapted to handling complex conditions (eg. non-homogeneities 
or anisotropy) as is the finite element approach, rt is noted that the 
finite element modeling approach can be used in precisely the same manner 
as the boundary integral method in that sub-regions are easily modeled 
provided the boundary conditions for the sub-region are available. Fi- 
nite element modeling is not restricted to sub-region analysis however, 
and is generally more flexible in its applications than other numerical 
techniques. 

General Description of Speckle Analysis 
With Finite Elements : 

Speckle interferometry is a technique for experimentally measuring 
"in plane" displacements (displacements normal to the direction of sight) 
of a loaded object. This technique makes use of the fact that the small 
imperfections of a surface uniquely "map out" that surface (the imper- 
fections determine a "fingerprint'bf the surface). To record the in plane 
displacements of an object, the characteristic pattern of the surface 
under consideration is obtained by illuminating the unloaded object with 
a coherent light source (such as a laser) and photographing the image 
of the surface. The object is then loaded, producing the displacements, 
and the photographic film plate is exposed a second time. The film 



now has a record of the characteristic pattern of the surface imperfec- 
tions and a record of how that pattern was distorted when the object 
was loaded. If a beam of coherent light is passed through a small por- 
tion of the processed film plate, the two patterns will cause the light 
waves to interfere and produce fringes which can be related to the in 
plane displacement of the illuminated point on the image. This experi- 
mental technique is well suited to the study of relatively flat members 
loaded in a plane stress state although other applications are possible. 

Once displacements are obtained, the conventional method of obtain- 
ing stresses is to numerically differentiate the displacements at the 
point of interest to determine strain components and then to use con- 
stitutive equations to obtain stresses. This requires displacements 
to be obtained at many points in the area of interest. Unfortunately, 
areas of interest (areas of high stresses - large displacement gradients) 
are not necessarily areas with large displacements. If displacements 
are below the resolution of the measurement technique, insufficient 
experimental data is available for local strain determination. 

Use of the finite element method allows (if necessary) displacement 
data to be obtained away from the area of interest and in an area where 
more accurate data can be collected. The finite element approach to be 
used with speckle interferometry should accept displacement data to spe- ■ 
cify the boundary conditions for the area under consideration. Hence, 
a conventional "displacement" formulation of the finite element model 
is suitable. The boundary conditions are determined by specifying the 
displacements surrounding the area of interest, not necessarily close 
to the points of interest. Additional boundary conditions, such as 
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stress free surfaces, are also readily accommodated. The finite element 
method then uses the physical characteristics of the structure along with 
the experimentally determined boundary conditions to determine the stress 
distribution throughout the area of interest. 

Description of Finite Element Program 
and Modeling Considerations 

A suitable element computer program into which experimentally ob- 
tained displacement data can be input has been developed in-house. The 
program employs a six-node, plane stress triangular element. The element 
nodal points are located at the verticies and midsides of the triangle. 

Two orthogonal components of displacement are permitted at each nodal 

point; thus, twelve displacements or, "degrees of freedom" are associated 
with each element. The displacement field within each element is interpo- 
lated from the nodal displacements through complete, second degree poly- 
nomial functions. Thus, linear stress distributions are permitted within 
each element. All criteria for guaranteeing convergence to the "exact" 
solution as the element mesh is refined have been met with this element. * 
The program is capable of modeling linear elastic, nonhomogeneous , aniso- 
tropic plane stress problems. Thermal stress modeling is also possible. 

Any combination of nodal point displacement values can be specified. Boun- 
dary tractions (normal and shear stresses) can also be specified as well 
as "point loads". 

To facilitate application of the program an automatic grid genera- 
tion subroutine was developed. By defining model boundary geometry an in- 
ternal mesh of triangular elements with associated nodal coordinates and 
connectivity parameters can be automatically generated with this 
routine. Figure 12 illustrates grid generation for a complex structure 
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Figure 12 


Automatic Grid Generation of a Complex Structure 
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by use of this program. 

This finite element program was validated on several trial prob- 
lems with known solutions and found to give excellent results. The linear 
strain triangle is relatively insensitive to aspect ratio and is far su- 
perior to a three node, constant strain triangular element. Rigid body 
movement can accompany the deformation without altering the. stress field. 
This feature is important since experimentally measured displacements 
usually include some contributions from rigid body movement. 

Some judgement is necessary regarding mesh refinement for problems 
where stress gradients are not known apriori. Consideration of the per- 
mitted stress variations across individual elements together with anti- 
cipated stress gradients in regions of stress concentration should form 
the primary guidelines for grid selection. 

When analyzing a structure or a sub-region of a structure by finite 
element analysis, boundary conditions on all defining boundaries of the 
grid must be specified before the problem can be solved. Either displace- 
ments (experimentally determined or otherwise), boundary stresses or com- 
binations of both can be used. The selection should be based on the bound- 
ary conditions that are known with the most certainty. If displacements 
are experimentally measured by laser speckle techniques, for example, it 
is likely that displacement boundary conditions can be determined along 
all boundaries. If a portion of the boundary is, however, free of applied 
traction, this stress boundary condition should be used in lieu of exper- 
imental displacement data along that boundary. This is because such a 
condition can be enforced and be totally free of experimental error. It 
is noted that most finite element programs are written such that boundaries 
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are automatically assumed to tie "stress free."' unless otherwise specie 
fied. Thus this type of stress boundary condition usually need not be 
explicitly specified with input data to the program. Displacements and 
non-zero boundary tractions are explicitly input as data. It should 
also be realized that random pertabations in boundary data (perhaps due 
to experimental error} can significantly affect local stress values. 
Saint Venant's principle indicates these effects to be localized in some 
sense provided that the boundary conditions are statically correct over- 
all . 

These comments are intended to provide some general guidelines for 
consideration when selecting boundary conditions from experimental data 

for use in finite element modeling. 

A copy of the finite element program is included in the appendix. 
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Application of Laser Speckle/Finite Element Technique : 

Two applications of the Laser Speckle/Finite Element Technique are 
presented. The first application is to a cantilevered beam with a con- 
stant distributed load applied from the middle to the end of the beam. 

The second application is to a "C-ring" in diametrical compression. 

The Cantilevered Beam: 

This example is given to illustrate how the completeness of descrip- 
tion of the displacements at the points of applied loads effect the accu- 
racy of the stress distribution obtained. Laser speckle data was not ob- 
tained. Instead, displacements were determined analytically from beam 
theory for an aluminum cantilevered beam 1 inch wide, 15 inches long and 
1/4 inch thick with a constant distributed load of 3 Ib/in applied from 
the middle to the end of the beam. For finite element analysis the beam 
was modeled using 20 plane stress elements, Fig. 13. 

Five different cases were analyzed. The first case [Case 1] had 
the displacement at the end of the beam due to the distributed load as 
the only input into the F.E. (finite element) program. The other points 
on the surface had their boundary conditions left to their defaul t val ues 
of being stress free (the stress free condition is the most accurately 
known boundary condition except in areas where a force is known to be 
applied, such as under the distributed load ),Case 2 had the displacements 
at the start and end of the distributed load input into the F.E. program; 
Case 3 had displacements input at the start, middle and end of the dis- 
tributed load; Case 4 had displacements input at six points (every other 
node) under the distributed load; and Case 5 had the displacements at 
every node (eleven points) under the distributed load input into the 
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finite element program. The normal stress, ax, at the surface versus 
the length of the beam has been plotted for these five cases and are shown 
in Figs. 14 - 18 . Note that as the displacements at the applied load 
were more completely described, a x approached the analytical curve for 
ax- 

For the second application it was desirable to apply the Laser 
Speckle/Finite Element Technique to a model in a nontrivial plane stress 
state. It was also desirable to chose a model for which an analytical 
solution was known. A "C-ri'ng" in diametrical compression satisfied 
both of these requirements. The stress state is nontrivial due to the 
curved beam nature of the model; an analytical solution is available in 
Theory of Elasticity , by Timoshenko and Goodier. 

The C-ring had an inside diameter of 3 inches and an outside diame- 
ter of 4 inches. It was made of aluminum and was 1/4" thick. 

Equipment was needed to hold the C-ring and to provide reference 
measurements. In order to load the ring in diametrical compression a 
special test rig was built, Fig. 19. This rig allowed displacement along 
a diameter of the ring by applying a force through a hydraulic cylinder 
arrangement. One of the main considerations in designing the test rig 
was to restrain the C-ring from any movement out of the plane of interest 
without inducing any additional forces which would effect the stress dis- 
tribution. 

In order to measure the diametrical displacement the C-ring exper- 
iences, a dial indicator (accurate to 0.0001 inches) was mounted along 
the diameter through which the force was applied. A calibration speci- 
men was also included ; n the experimental setup so that an object under- 
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*** Quality 


Figure 19. Load Arrangement For C-Ring 



going a known rigid-^body displacement would be recorded on the speckle- 
gram. 

Speckle interferometry was used to determine the displacements which 
the C-ring underwent during the diametrical compression. The model was 
illuminated by a helium-neon laser beam diffused through a spatial filter. 
The images of the model were recorded on a Agfa-Gevaert 10E75 4" x 5" 
film plate held in a camera arrangement. The first exposure of the film 
plate recorded the image of the ring having a small amount of preload 
applied to it. The second exposure recorded the image of the ring after 
applying a diametrical compression which reduced the diameter along the 
line of action by approximately 0.01 inches (referenced to the preloaded 
position) . 

Data was taken from the speck! egram along two radial lines (A-B 
and C-D on Fig. 20) separated by 30° on the C-ring. Spacing and orienta- 
tion of the Young's fringes were measured along each cross section at 
13 points correspond!' ng to the nodal points on the finite element model. 

Young's fringe information was recovered from the speck! egram by 
illuminating the film plate with an argon laser whose beam was converged 
with a long focal length Tense. The halo was projected onto a scribed 
piece of frosted glass mounted such that it would be rotated to allow 
alignment with the fringes so that their angular orientation could be 
measured. 

Accurate placement of the speck! egram in the laser beam was accom- 
plished by using a translating platform which could be positioned in 
three directions with micrometers. This platform also allowed rotation 
about the axis parallel to the laser beam for ease of determining a 
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reference orientation for the Young's fringes. 

Due to symmetrical considerations, only 90 degrees of the C-ring 
(from the point of application of the load to the center) needed to be 
modeled for finite element analysis. This portion of the ring was divid- 
ed into 216 elements with the density of the elements increasing toward 
the inside edge, (Fig. 21). The increased density of elements close to 
the inside edge was necessary due to greater stress gradients in this 
area. 

To ascertain that this finite element grid would yield accurate 
results, the principal stresses from the finite element program were 
compared to those from the analytical solution. These stresses from both 
solutions along the cross section in the middle of the model are shown in 
Fig 22 . The largest difference between critical stresses of both solu- 
tions was less than 1.02» using this grid. 

The Young's fringe data taken along the 0° and 30° lines were re- 
duced to obtain displacements in the x and y directions,. These displace- 
ments were plotted and a "best fit" straight line was drawn through the 
data in order to smooth experimental error effects, (Fig. 23 and 24). 

Data taken off of the lines were used as displacement inputs in the fi- 
nite element program. These displacements defined the boundary condi- 
tions along the 0° and 30° lines; other boundaries (inside and outside 
edges) were specified as being stress free. 

The validity of using a "best fit" straight line for the displace- 
ment data was investigated by plotting the analytic displacements along 
the 0° and 30° lines. A straight line drawn through this data showed 
very little deviation of the data from the line. Using the displace- 
ments taken from these lines as inputs, the finite element program 




Figure 22. Principal Stress Along 15° Cross Section Of C-Ring 
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determined stresses which were in very close agreement with the analy- 
tically determined stresses at points removed from the cross-sections 
where the displacements were input. 

Principal stress results along the 15° radial line from the finite 
element program using the straight line fit of the experimentally deter- 
mined displacements is shown in Fig. 25 . It should be noted that the 
stresses have been scaled such that the Speckle/F.E. determined princi- 
pal stress on the inside edge match the analytically determined princi- 
pal stress at this point. This scaling was done because an experimental 
scaling factor could not be determined due to the fact that the calibra- 
tion specimen appeared to have undergone other than a purely rigid body 
in plane displacement. This only limits comparisons of the magnitudes 
of the stresses and does not hinder the comparison of the trends of the 
stresses . 

Examination of Fig. 25 shows that the analytic and Speckle/F.E. 
principal stresses varied from the forced difference of 0 psi on the in- 
side edge to a difference of approximately 1400 psi at the outside edge. 
Note that the shapes of the curves are very similar and only their place- 
ment and orientation vary significantly. 

Figures 26 and 27 show new straight line fits through the experi- 
mental displacement data. Although these new straight line fits do not 
seem to be "best" fits, they do appear to be within the range of un- 
certainty for the experimental data as is indicated by the scatter of 
the data. When displacements from these new lines were input into the 
finite element program there was much better agreement between the prin- 
cipal stress curves as can be seen in Figure 28 . The difference between 
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Figure 25. Principal Stress Along 15° Cross Section Of C-Ring 
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the Speck! e/F . E - and the analytic principal stresses varies in this 

case from the forced 0 psi at the inside edge to approximately 350 psi 

* 

at the middle of the cross section. Such a large variance in results 
for the small change in displacements input dramatizes the sensitivity 
of the stress results to the experimentally determined displacements. 

The large scatter of experimental data about a straight line (as 
opposed to the straight line nature of the analytic data) indicates high 
uncertainty in experimental data. The fact that the Speckle/F.E. method 
was able to determine the general trend of the stresses from this un- 
certain data shows definite promise for this method. A simple numeric 
differentiation of displacement data would! have yielded stresses not 
very representative of the actual stress state. However, in cases simi- 
lar to the "C-Hng*" accurate experimental determination of displacements 
within a straight line approximation should allow the Speckle/F.E. method 
to closely determine the stress state. 
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Appendix 

Finite Element Program 


r * 

r. 'MNSTSGO rs A fenitf element STPJCT'JPAL analysis PROGRAM written * 
r BY IP, .JOHN l t T'JRN r R , ASSOCIATE PROFESSOR OF AF-R I CULTURAL ENGINEER- * 
TNG AT AIJPtJPN UNIVERSITY. PLNSTSGD USFS SIX NGOFD PLAIN STRESS TR I - * 

ANG fl AR FLFMFNTS. FCNTAINEO IN PLNSTSGD IS A GRID ROUTINE? DEVELOPED * 
F OY FRANK H, VALADE, JR. AT AUBURN UNIVERSITY, WHICH GENERATES FLE- * 
r. MFNTS FOR A PASTC GEOMETRY OFSCRIBED BY THE USEP . * 

<■ * 

r 


CGGGRGGGGGGGsnGGSGGGSGGGGGGGCGSGGGGGGSGGGGGGGGGGGGGGGGGGSGGGGGGGGGGGGGGG 
r A/14/SX THIS PROGRAM HAS BEEN MODIFIED FROM THE VERSION IN STRSLI3. 
r THP FHANGFS ARE SET OFF BY ROWS OF AM»EF SANDS ( G ) . JOHN WEATHERS 
C£££FG r . ££££££££ £££££££££££££££££££££££££££££££££££££££££££££££££££££££££ 


IMPLICIT RFAL«-3(A-H, C-Z ) 

DIMENSION XG (4*1) »YG( ABU , IMATL (21 M . 0 ( 1 , 3 , G ) , NO C ( 2 1 6 , 6 ) , B ( 962 , 
*54) ^lo^l 

0 I MENS TON ST D NX ( 216,6), STRN Y ( 21 6 » 6 ) , STRNX Y( 216 ,6 ) 

DIMENSION STRSX (216 ,6) ,STRSY{ 216,6 i ,STPSXY{?). 6,6 ), TITLE! 10 ) 

NR 0 = 5 


NWT = 6 

^ £££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££ 
READ {NR 0.5 00 ) T I Tl F , TH I CK ,R Bt) I SP 


500 FQRM AT ( l 0A° , /,F1G.0,F1 C.O) 

WRITF(NWT, 600) TITLE. THICK, RBD ISP 

AOO FORMAT { • l* » 1 0 A G , / » • 0 * , * TH ICKNES$= ’.FlO.A./.'C’t’R.B. DlSP=',F10.4) 


F F.££££F £££££££ r.££T ££££££££ GS £££££££££ G £££££££&£££££££££££££££££££ 8 ££££££ 


RFAOINRO, 1 )NNP,NFLE,NMP t NRnw,NCCL 

l rf)R M AT (5 15 ) 


V.R I TE ( NWT , ?) N N P » N F L F ,NM P * NROW , NCCL 
7 t^ORMATf i O' , 'NNR=« , T 5, /, 'C' , ’NELE= ' ,1 5, 

* / , * 0* , »NMP= • , I 5, / , ! O' , »NRPW= » , 15, / , • C* NCOL = ' ,15) 
CALL 4MA I N ( NP D, NK T , NN p , N ELE , NMP , NRCW , NCCL , X G , YG , I« AT L, 0 , NOD , 

* ?,F,STRNX,ST p -NY, STPNXY, STRSX, STRSY.STRSXY. T I TLE, THICK, R 801 ) 

CALL EXIT 

STOP 


FND 


, 
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SUBROUTINE AMA !N( NPO, NWT , NNP , NELE , NMP ,NROW, NCOL * XG , YG, I MATL , Q , NOO , 
*3, F, STRNX, S T PNv t STRNXY, STPSX, STRSY , STRSXY , TITLE. THICK, RBO ISP ) 

r 

IMPLICIT PEAL*3( A-H,0-Z) 

DIMENSION XG(NNP) ,YG(NNP ) » I MATH NELE) , C< NMP , 3 ,3 ) , NOD (NELE , 6 ) , E ( 3 , 3 
*) ,S( 12.1 2) , 3(NROW,NCOL ), F(NROW) 

DIMENSION M3CL (100) , VO ISP (100), STRNX (NELE, 6) , STRNY ( N EL E ,6 ) , STRNXY 
* ( NELE, 6 ) , STRSX(NELF,6) , S TR SY ( NE LE , 6 ) , STRSXY I NELE ,6 ) , TI TLE ( 10 ) 

LONG =100 

CALL GRID01 ( NNP , XG , YG , NEL E , NOD ) 

WPTTEfNWT, 17) 

17 FDPM AT (•].', 12X, 'NOOF N JM BER ' , LA X , ' XG ' , 2B X , * YG ' ) 

00 7 \=i f NNP 

** 0 E AO ( NRO * 1 ) I ,XS( l ) ,YG( t ) 

7 WRITFKNWT, 1.6) I , XG ( I ) , YG ( I ) 

16 FORM AT ( ' O' ,13X, I 5 , 1 5X , F l 0 .4 , 20X , F 10 . A) 

1 c 0 R M AT (1 S t ?F10.0) 

0 GGGRGSSEESSGSS&S5GSG&&S&&G&GG&GEG&G&GaGGGGGG&GGG&E&6&&G&&&&G&G&&& 

C P E AD ( N R D, 7 ) ( IMATL ( I ), 1 = 1, NELE) 

on 100 I=L»NFLE 
100 IMATLU1-1 
C WRITE(NWT, 1 9 ) 

IQ cqrmaT ( • L* , 1.0X) 

r WPITF(NWT f ) 3 ) ( I,IMATL(I),I = 1,NELE) 

IP FORMAT! *0*', ' l MATH ',13,' i=' ,15) 
r GGGGGGGGGGGGGGGGGGGGGGGGGGG&G6GG6GGGGGGGG&G&&GGGGGGGGGGG&6GGGGGGG 
? Ff)RMAT( 5 0 1 1 ) 

n Q 1 I=l,MMO 

B P EAD(NRD,A)Q( 1, 1 ,1 ) ,C ( I ,1 ,2 ) ,Q(I »l ,3) ,Q( l ,2 ,2 ),Q( 1,2,3) ,0( 1,3,3) 
WRT TE(NWT t ) 9 ) 

WPITFINWT, 29) (((I,J,K»Q(I»J,K) ,K = J,3) , J= l , 3 ) , I =1 , NMP ) 

?0 FCRMATC O' , • 0 ( • , ll , • , ' , 11 ,• , • , II , » ) = * , E14„7) 

4 F C RM A T ( 6 E 1 0 , 0 ) 

^ nn 5 N = 1 , N p L F 

r 6 R EAT ( NR D , 6 ) I , ( NOD ( I , J ) ,J = i»6) 

6 FORMAT (7 15) 

WRITF(NWT,7l ) 

?1 FORM AT { ' 1 ’ , * FLEMENT N UMB ER » , 20X , » NOD E NUM3FR S' ) 

DO 1=1, NELE 

30 VsR ITE(NWT ,?? 1 I , <NCO( I , J) ,J = i ,6) 

27 FORM AT ( »0»,5X,T5,10X,6(I5,5X)) 

On io 1= 1 , NP QW 
^9 10 .1 = 1 ,NCCL 
in B ( I , J )=0 .0000 
on r I=l t NFLF 
N=IMATL( I ) 

CALL MATL (N,Q,NMP,F) 

I F ( N . EO • 1 8 ) W D ITF(NfeT,13) E 
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CALL STTFF( XG, YG,NOO,NFLE»NNP,S»E, I .THICK) 

CALL ASS CMP ( B.NROW.NCOL, S.NOD.NELE, I ,12,6,2) 

WRITFINWT.6) I 
R CONTINUE 

WPITf=(NWT,n) ( (B< I, J) , J = 1 »NCOL ) ,I = l,NROW) 

WR I T E ( NWT ,601 M F ( I ) , 1 = 1 ,NROW) 

CALL BCPnRC(NRD,F f NRQW ) 

WP ITF(NWT, 500) ( F( I ), 1= l ,NRO W) 

500 FORMAT* /,] OX.6F15.4) 

CALL RCSTF3(R, c ,NP0W,NC0L ,MBCL, VD ISP ,LCNG , NC I SP , NCOL yRBOISP) 

WR I T E ( MW T ,1 3 ) { { R ( I , J ) ,J = 1,NCGL) , I = l,fsRCW) 

WRITE* NWT, 601 ) (F (I ) , T=L ,NRDW) 

CALL SYR AT ( B,F, NROW, NCOL , 1 , NRG W ,N CO L ) 
cQRMAT(// f ( / , 6 F?0 • 6 ) ) 

WRITE (NWT,P5) 

75 FORMAT ( M • , RX, * NODE ' , 14X, * X-DI SP. • ,22X, • Y-OI SP. • ) 

OQ 1 = 1 ,NN|P 

0 = 2* I 
K = J- 1. 

WR ITEtMWT, ?4) I,FCK) ,F( J) 

P4 FORMAT (* O' , 15, 15X, F14.7, 15X.E14.7 ) 

RR CONTINUE 

WRITE* NWT, 7 7 ) 

? 7 FCRMAT ( • l • , • FLE . NO.*, 3X, 'NODE NO. * , 8X , * S IGMA-X * , 1CX , * SI GMA-Y • , 
*1QX, * TAU X- Y ' , 10X , • EPS ILCN-X' ,1 OX , ' EPS ILCN-Y • , 10X,» GAMMA X-Y* ) 
rjn i/ + T =lfNP |_p 
N ~ I M AT L ( I ) 


CALL MATL(N,0,N«PtE) 

CALL STP AIN (XG,YG.NOO,NELE,NNP,E, I , F , NRG W » S T R S X , STRSY, 

*ST°.S XY, STRNX, STRNY, STRNXY) 

00 15 K= 1 ,6 

WRITE (NWT.50P } [ . NQO ( I , K ) , ST RS X ( I , K ) , STP S Y ( I , K ) , S TR SX Y { I , K ) , 
♦STRNX <T , K) , STRNY ( I ,K) .STRNXY ( T ,K) 

50P FORM ATI 'O’ , PX» T?, 6X, J4, 1 1 X . 6( E14.7, 3 X ) ) 


C £ GGGGG ££&££££ GGGGTf G ESC GGG CSC GGGGGGG GGGGS GGGGGGGG GG GG GGG G GGGGGGGG 


STSMM = OSORT ( { (STPSX ( I ,K ) -STRSY ( I ,K ) ) / 2 3 0 000 ) **?. +ST RS XY ( I ,K )**2 ) 
STSMM1 = 0.5D00*( STRSXIl ,K) +STRSYM ,K) J+STSMV 

STSMM? = 0. 5000*1 STRSXI I ,K) +STRSYI I ,K) )-STSMR 
W R IT c ( NWT , 5 OP ) STSMM] ,STSMM2 


f OR FORM AT(] X,T?‘7, P( E14 .7 ,?X ) ) 


FGGGGSG GGGGGGG GGGGGGGG SCGGG GGG GGGGGGGG GGGGG&GGG&&G&GGGG&GGG&&&GGGG 


if rrNT TN'JF 


14 CCNTfN'UF 


R F tijr M 
CNC 


1 


► 
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r, 

r, THIS ROUTIN' - * ASSEMBLES A SYMMETRIC BANCEC STRUCTURAL STIFFNESS 
" MATRIX VIA INCtVfnuAL FLFME N T S . 


SUBROUTINE ASS FMR ( R » NRQ ,NCO , S »NOD , KFLF , NTR I . NO I M . NNP » NOFPN) 

C 

r ****** ******* ************* ********************************************* 

* 
* 
* 
* 
* 
* 
* 
* 
* 
# 
❖ 
* 
* 
* 

^ ************************************************************* ********** 


C 

»■> 

r 

r 

r 

r 

#* 

* 

r 

c, 

r 


P * GLOGAL STIFFNPSS MATRIX! SYMMETRIC BANDED) 

S = ELFMFNT STITPNFSS matrix 

NNP = MP CF NTH At POINTS IN ELEMENT 

N DF D N = NO OF DEGREES OF FREEDOM REP NCDE 

NOIM = D I MFNS TON OF FLEMENT STIFFNFSS MATRIX = NN p *NDFPN 

NOD = GLOBAL CCNNICTIVITY VECTOR 

NTR I = p L P M P N T NO 

NRO = NO OP ROWS IN STRUCTURAL STIFFNFSS MATRIX 

MCOL = NO OF COLUMNS IN TOTAL STRUCTURAL STIFFNESS MATRIX 


r 

IMPLICIT RFAL AC HA-H.C-£) 

DIMENSION NQC ( N FL P t 6 1 .B ( NRC *NCC ) * S tNCIM.ND IK) 
c I , T N index ROWS I, UN TNOFX COLUMNS 

On 1 I si, NN» 

NR ° W = ( NO n ( NTR 1,1) - l)*NDFPN * 1 
TO I J-1 .NNP 

Nrni a ( NOD! NTR l » J I - l > *NOFpN +• l 
r CHECK Tn S n P IF RFLCW MAIN DIAGONAL 

IF! NROW. GT.NCOL) (.0 TO 1 
MBCOL a TARS! NCOl-NRCW ) *1 
r OO PN NUMBER TF OFGREES OF FREEnOM PEP NQOF 

r>C P I N« I .NOFPN 
OP 3 s JN=l. NOFPN 

C. CHFCK POR EL C *1ENT S BELOW MAIN DIAGONAL*. 

I F (NBCOL . EO. 1 ) GO TC 10 
GO to H 
10 CONTINUE 

IF (JN .LT .IN) GO TO ? 

U CONTINUE 

NRCPll = NBCOLfJN-TN 
N ROW l a NRPW-L*- IN 
MRS = N n p PN« ( I -I ) + IN 
NCS a NPFPNAf J- l ) + JN 

R (NROWl .NBCOl. 1 ) = P ( NROW l . N BCOL U F S(NRS.NCS) 

*> continuf 

1 CONTINUF 
RFTUPN 
FND 
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SUBROUTINE RCFORCI NRB,F, NROW) 

IMPLICIT REAL*B( 4-H, 0-Z) 
niMPNStON F(NRQW) 

00 2 T=1 ,NROW 

F ( 1 5=0.0000 

READ (NRC, UNFORCF 

VP IT F ( 6, M ) NFORCE 

FORMAT! • I * , IQX » ' NFORC E- ? ,15) 

JFINF0PCF .FQ. 0) RETURN 

R EAD ( NR 0 1 1 1 INVAR, F INVAR ) ,N = I, NFORCE) 

FORMAT I 51 IF,E10.4> ) 

no 10 1=1, NRGW 

IFI FI II . EO.O > 00 TO IC 

WRITFI6, A) I , F ( I ) 

FORMAT! * O' . »NVAR=« , 15 ,10X,' FORCE=* ,E14 a 71 

CONTINUE 

RETURN 

FNO 
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S JBPO JTINF BCSTF 3 IB, F f NV,NC,MBCL,VDISP,LCNG,NOISP,NeAND f RBOISP) 

r 

I M D L TC I T RFU* 9 < A-H.C-Z) 

niMPNSION R ( NV » NC I » F ( NV» 1 ) • MBCL ( LONG ) ,VD ISP (LONG) 

C 

C ** **** *££*#** lk#*#)}!#***#****#* * ***# ******* $ * $ * + + *j(t Jjc$ $ $ # j#ea>: $ sfr# $ 

r 

C THIS ROUTINE PLACES NCN 7 .ERC C I S P LAC EMENTS INTO THT SYMMETRIC STIFFNRS 

C 

r *********************************************************************** 

NRD =5 

p HAH ( 5 , 100 ) MO ISP 

1 00 FORM AT (41 10) 

ND.TSP = MDISP 

r. NO t S D = NO OP S p C C I F T EO DISPLACEMENTS 
C T HFS p INCLUDE PIGIO BOPY DISPLACEMENTS 

p FAD ( 5 , 101 ) (VDTSP(K) ,MBCL(K) *K= 1 ,M 0 ISP) 
CEEEGEEGEEEEEEEGESEGEESSEEEEGEEGESGEEEEGGGEE ECefiKCSS 6 £GS 66 E££fiefi 6 SS£Efi 

101 P OR M AT(F] 0 8 4 , n 0) 

101 F0PMAT(4(Et0.7,I5)) 

C 

r ***** P« IT IN RIGID BODY DISPLACEMENT ***** 

C 

no looi k = i » mo i s° 

1001 VOTSP(K) = VDTS°(K) +-R 3 DI SP 

r 


C 6 SEE EG EGG EG EG GGGGGEGF, CGEEESE GEEEEESSESEEEEGSGGEGEEEEEEEEEGEEEEEGEEEEE 
C MRCL = NO OF VARIABLE AJ WHICH DISPLACEMENT IS SPECIFIED 
r v n T c, P = VALIJP OF DfS PLACEMENT 
WR I T F ( 6 » 10?) WO ISP 
in?. FORMAT! ’ 1 * , 1.0X, •NOISP=', 15) 

WRITE (6, LOOM MBC L ( K ) , VDI SP(K) ,K=1,NDISP) 

1 03 FCRM AT( • O' , 'NVARt* t 15, 10X *• DISPL. = ' * El.4.7 ) 

DO 1 1=1 , NO I S D 
N V AP =M BC L ( I ) 

VAL=vni$o{ i i 


00 2 J = 1 , MB AND 
NP=mv AP *■ l- J 
NRM = NVAP+J-1 
r F (NR — 0) 3.3,4 

4 CONTINUE 

F(NR,1) = F ( NR , l ) -VAL*B(NR,J) 

R (NP ,J? = 0.0000 
n CCNTIMUE 

IF(NPM-NV) 5,5,6 

5 CONTINUE 

F ( NR M , 1 1 = F ( NRM , 1 1 - V A L*B { NVAR , J ) 

6 CONTINUE 
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B(NVA°,J) = C . 0000 
7 rnNTIN'Jr 

R ( MVAR « l ) = 1.0DC0 
R ( N VAR » L ) = VAL 
1 COMTIN'JP 

PFTUP.M 
FNO 
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!! 
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f 
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SUBROUTINE OMINVI A,N,D,L*M) 


• • • ••••• •••«•»■ • ••••••••»• • • 1 * •••••••••. ft o • • • • • » • 

SUBROUTINE OKINV 
PURPO SF 

INVERT A MATRIX { OOU3LE PRECISION VERSION) 

I SAGE 

CALL DMINVI A,N,D,L,M) 

OFSC.R IPT ION CF PARAMETERS 

A - INPUT MATRIX, DESTROYED IN COMPUTATION AND REPLACED BY 

PFSULTANT INVERSE. I 

N - ORDER CO MATRIX A | 

D - RESULTANT DETERMINANT i 

L - WORK VECTOR OF LENGTH N j 

M - WORK VECTOR CF LENGTH N j 

/ 

REMAP KS j 

MATRIX A M l ST BE A GENERAL MATRIX 

[! 

r SUBROUTINES AND FUNCTION SUBPROGRAMS REOUIRFO 

C NON- !| 

r || 

0. M FT HO D jj 

0. THF STANDARD GAUSS-JORDAN METHOD IS USED. THE DETERMINANT 

C T S ALSO CALCULATED. A DETERMINANT OF ZERO INDICATES THAT 

C THP MATRIX IS SINGULAR, 

r 


r 

IMPLICIT p fal*B( a-h, c-z ) 

O IMF NS ION Ml) ,L(l)fM( 1) 

r 

C SEARCH FOR LARGFST ELEMENT 

C 

0 = 1.0 DOO 
NK=-N 

n n BO K= 1 , N 
MK = MK 4-N 
L ( K )=K 
m (K) =K 
KK=NK+K 
BI G A = A ( K K ) 

DO ?0 ,J = X,N 
t Z=N*( J- l ) 


n n «"> o 
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r 

r 

r 


no ?o i=k,n 

T J= I 7 + I 

10 TF{QA8S(BI0A)-DA8S{A( IJ ) ) ) 15,20,20 
15 3IGA = A ( I J ) 

L ( K ) = I 
M (K ) = J 
20 CONTINUE 

INTERCHANGE ROWS 

J =L ( K ) 

I F ( J -K, 1 25,35,95 
25 KI=K-N 

on 2 0 1=1, N 
Kl=Kl+N 
HCLO=-A( K I 1 
J I = K T-K+ J 
A ( K T ) = A ( JT ) 

2 0 A ( J I ) =HPl_n 


0 INTERCHANGE COLUMNS 

r 

25 I = M ( K ) 

IF(I-K) 45,45,38 
28 JP=M*(I-1) 

oc 40 J= 1 , N 
JK = NK f J 
J l=JP-*-J 
HCLO=-A( JK ) 

A ( JK) = A( JI ) 

40 A ( J l ) = HCL 0 

DIVIDE COLUMN BY MINUS PIVOT (VALUE OF PIVOT ELEMENT IS 
FONT A INFO JN 8 l c A ) 


45 IF(RICA) 48,45,48 
A A 0=0.0000 
RETURN 

4« on 55 1 = 1 , N 

I F ( I — K V 50,55,50 
50 T K = NK+ I 

A ( T K ) —A ( IK)/ (-2 IGA) 
55 CONTINUE 

RFOJCE M'TRIX 


I 


on 65 1=1, N 

IK=NK+I 


U C 


65 


i 


f 


C 


r 

C 

C 

r 

r 


r 

r* 
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HrLO=A( IK ) 

I J= I— N 
1C 6 5 J = 1 ,N 
T J=T J4-N 

IF(l-K) 60,65,60 
60 IF(J-K) 6?, 65, 62 
6? KJ=TJ-H-K 

A(IJ)sHnn*A(KJ)tA( IJ1 
65 CENT INIJE 

OIVtOF ROW RV PIVOT 

K J = K-N 
on 75 J = l » N 
K J=K J+N 

tF(J-K) 70,75,70 
70 A (K J ) =A( KJ1/BIC-A 
76 CONTINUE 

PPOO'JCT OF o IVOTS 
0 =D*BIGA 

P E°LA C E DIVOT BY RECIPROCAL 

A {KK 1=1. 0000 /RIGA 
«0 CONTINUE 

FINAL ROW A NO COLUMN INTERCHANGE 


K =N 

100 K=( K-l ) 

T F ( K ) 1.50, I5C, 105 
105 I =L ( K ) 

T F { I — K ) 120,120, IC8 
100 JQ=N*lK-l) 

JR=N*( t-1) 

00 110 J = l , N 
J K = J 0 J 
HCLO= A (,JK ) 

J I -JR + J 
A ( J K 1 =— A ( J I ) 

110 A ( J T ) =HOL0 
i .?0 J=M(K) 

IF(J-K) 100,100,125 
175 K T = K-N 

DO 130 1=1, N 
KI=K I*N 
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HCLO=A(K 1 ) 

J I=K I-K+J 
MKt) =~A ( J T ) 


130 4(J!) =HPLP 


00 TO 100 

ORIGINAL PAGE IS 

150 RETURN 
ENn 

OF POOR QUALITY 
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t 


t 


SUBROUTINE GPTOOU tl.XO.YG. L2,NCD) 

IMPLICIT R FAl *3( A-H»n-Z) 

DIMENSION p $PAN(?0),NnD(L2 f 6),XG(Ll), YG(L1 ) 
PEAC( S, 15) N p AP 
READ(5,15) NUMS 
WRTTF(6.20)NO*R,NUMS 
°0 FORMAT ( *!*,?( ) 

1 5 FCRM AT ( 1,615) 

CO l T= It NUMS 

1 PFAR(5, 16) P S' 5 AN ( I ) 

WR I T F ( 6 t 7 1 MPSPANU) *1*1, MUMS) 

21 FORMATt *0 • ,5( F10.A, 2X) ) 

K= 1 

1 1 = 2 #N' J MS 4-2 
DO 2 J=1,MPAP 

RFAT(5t 16) STRXt STRY, STPXtSTPY 
WRITFf6t7l) STPXtSTPY.STPXtSTPY 
FORMAT ( 8F10. A) 

XG(K)=STRX 

Y 0 ( K ) = S T R Y 

SP ANX=ST° X-STR X 

S°A N Y=S TPY-STR Y 

CO 0 L=ltNUMS 

FR AC=PS P A N ( L ) / 100. CHOC 

FNC X=SP ANX*FRAC 

F NC Y = S P A N Y * F R A f 

K i = K+ 1 

X0(«) )=X0< K )+FNCX/2,0000 
Y G ( K 1 ) = YG { K ) f C NCY /2 • OOC 0 
K 2 = K + 2 

XG ( K 2 ) = X G( K ) 4-FNCX 

Y G ( K 2 ) = Y • G ( K) f F N C Y 
2 K =K 2 

2 K = K + II 
K = I l 

N0 = IT- l 
Nl=MDA p -I 
DC 4 J= 1 ,M1 
k°=k+no 

KM=K-NO 

XG(K)=0.5000* ( XG( KM ) f XG (K° ) ) 
YG(K)=0,5nOC*(YG(KW)+YG(KP) ) 

CO 5 L=l,NiJMS 
K l=K + 1 
K.w = KM + 2 

XG( Kl) = 0.50 00* ( XG< KM H-XG(KP ) ) 

VGfKl )= 0.50004 ( YGf KM )+YG(KP> ) 



8 

4 


7 

6 


r q 

c 

0 Q 
C 17 
1 3 
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K2=K+7 

KP-KP+? 

XG(K?) = C.SDOO*( XG<KM > + XG!Kp) 1 
YG( K2}= 0.5000*! YG< KY ) +YG1KP) ) 

K = K 2 
K=K+I1 
10 = 0 
1“ L 

12=11-2 
CO 6 K=l,Nl 
I A= T C*N0 + l 
18= ( rc+i)*No+i 

If = ( 10 + 7 ) >fefvJC + 1 

no 7 j=i ,nums 
NO 0 ! I , l )= IA 
NOO { 1 , 7 ) = I B 
None 1,81 = 10 
NOO ( I ,41=10*1 
N00< 1,0 1 = I A +.? 

NOO ( 1 ,h >=IA+1 
1 = 1 + 1 

NOO ( 1,1 > = I C 
NOO { 1,2 )= IC + l 
NOO ( I ,8} =IC+7 
NOC( 1,4 1=18+7 
NOO ( l , 0 1= [ A+? 

NCO ( 1,61=18+1 
IA=TA+? 

18=18+7 
IC = IC +7 
1 = 1 + 1 
10 = 10 + 7 

WR I T c RESULTS 
L l = ( N°A P *2- 1 1* (2*MJMS + 1 ) 

00 8 N= 1 , Ll 

WRITE(6,17)N,XG(N1,YG(N) 
L2=?*NUMS*(N D AP-1 1 
on o N = l ♦ L 7 

WRITE! 6,1 8 )(N, ( N GO ( N , J ) , J=t,6) ) 

FORMAT! »0*,I5,2(8X,F10.4)) 

FORVAK iq' ,715 1 

RFTIJRN 

FNO 
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5. UPROUT TK": M ATL ( f\ , Q , NM° f g ) 

IMPLICIT A-H f n-7J 

HIMFN^tnN QCNMn,^) ,01 <3*31 .F(T,T) 

•T l !l*l,^ 

TH 1. 

ni( 1 1, Jl lalMN, Il.JII 
1 T l ( J 1 » 1 1 ) =0 1 { 1 1 1 J t I 

i t - 0 ui,nMoi(?t?i*ciu. r n-in(2,'*)**?)-iu<i*?)*(Qm»at*Qin t *u- 
*Ql( ?,1)«Ql( 1 ,’'11+011 1 , ( l» 2 )*Qi(?,?)-Ql (2 ,3 )*Ql (l,*) ) 

i* ( i ♦ i ) it p , ? ) *oi n, ? ) -o u ? , t ) **? 

f (1,2)=-. (31 (l, •>1*01 n^)fOL( l,3)*QU2,^) 

r I ^ ^ l 1 ( 1 » 1 ) *Q 1 ( ? , 5 )-Q l ( l » ? )** 2 
ic * JIM,* 
m i Kl«* Jl, * 
r (Jl,KM«F< Jl,Kl)/OT 
r (Kl , Jt )»*=( JI ,Kll 
1 (TNT T NUT 
pet-irn 
r Nn 
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s JBP.nijT I N F ST IFF( XG, YG,NOO,NFLE ,NNP , E , I E , TH ICK ) 

TMPUCIT »FAl *S( A-H,C-Z) 

1 tMFMS ION XG( N\J° ), YG(NNP) ,N0D(NFLF,M,S(12,12), 

^■r(3»^)»X(6) »Y(6),C(6,6) 
p IMF NS KIN XX I NT < 6, 6 ),Y VINT (6, 6) , XYINTI 6,6) 
niMFNSIPN m(4),TT2(6) 

KNt=NQO( IF, l ' 

KN?=NOD( IF » 2 ) 

KNT = N00( IF,?) 

KN4=NOn( IF, 4 ) 

KN*5 = NQD( TF, 5 ) 

KN6=N00{ IF.M 
X (U = XG{ KN1 ) 

X ( 2) =XG( KN? ) 

X {? ) = XG( KN"* ) 

X (4) = XG ( KN4 ) 

X ( 5) =XG( KNS ) 

X(M=XG( KNM 
V ( l ) — YT < K N l ) 

Y ( 2 ) = YGf KN? ) 

Y ( 3 ) = Y G ( KN^ ) 

Y (4 )=YG( KM4 ) 

V ( 5) =YG( KNF ) 

Y (M=YG( KN6 ) 

TO l 1=1 ,6 

c ( i , n=i . oooo 
C ( I , -> ) =x f I ) 
r ( i , ?)=v ( i ) 
r. ( t,4)=x ( n yy(D 
r < t ,5)*x m*xu ) 

1 C( l ,S1=Y( I)*Y( I ) 

CALL f)MINV(C,ft,n,m , IT2 ) 

XC-(XU) *xm«-Xm)/l.QDOO 
YC=(Y(l) +Y( ?)>Y(«n )/^.OOCO 

ro ? i-i ,6 

X ( I )=X(! )-XC 
*> YIIUYU )-YC 

A As 0. 5000* (X(l)#(Y(3)-Y(S)M-X(?)A(V(S)-Y(l)1 + 

*X (F) *(Yt T ) — Y ( 1) ) ) 

A = C A B S < A A 1 

p l=XC*Yr +( 1 .0000/ 12.CDC0 )*(X(1)*Y(U*X(3)«Y(?)<-X(5)«Y(5)) 

D 2=YC#YCM l. 00 00/ 17.0000 )*( Y< l ) * Y ( \ ) «-Y ( 3 ) *Y ( ^ H- Y ( 5 ) * Y C 5 ) ) 
PT=xr*xr +n .oooo/ip.onoo )*(x(L)*xm+x(3)*xn) + x( , >)*x(6n 

PO 3 1 = 1,6 
OF 3 J=l ,6 

XX INK T, J )=A*{C( 2, I )*C<2, J) f2.0D00*XC*(C< 2, I)*C< 5,J)*-C< S, I )* 
*f ( ?, J ) )+YC*(C { 4, l . )*C< 2, J )+C< ?, I ) *C( 4 ,J) )+2.CnOO*P1*( C( 5 ,1) * 
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*ci4, j h-ci a, n*c is , j ) h-ci a. n*ci4, j >*P2+4.oooo*c( 5, n*c<5, j >*P3 ) 

3 YYJNTU, J ) = A* (C( n*C(3t J )«-2.0DC0*YC*(C( 3* l)*C( 6 , J ) *C l 6,n* 

*F (3, J ) )4-XC*( Cl 3, I)*C.( 4, J )+C<4, 1 ) *C( 3, J) ) «-2.0n00*Pl*(C(4, 1) *C(6 , J) ♦ 
*C(6, I)*C< A, J) ) *C (4, n *C ( 4, J »*P3+4.0D00*C (6, ! 1*C ( 6, J )*P2 ) 
no 4 1=1,5 
K = I + 1 

30 4 J=K , 6 

XX [NT (J, I)=XXINT( I, J) 

4 YYINT (J, I ) = YYINT( T , J > 

on 5 1=1,6 
on 5 j=i,6 

XYTNTU , J) = A*(C (2, I 1 *C ( 3 , J ) +XC* ( C ( ? , I ) *C < 4, J. ) +2 , ODOQ*C < 5 , I )* 
*r<3tJn*YC*(2.3n00*C(2,I)*C<6 t JMC(4,I)*C(3,J)>MC{4fl)*C(4,J)«- 

*r (5, n*c( 6f J m.onoo 1 ^ 1 + 2 , qdgo*ci 4, ti *c< 6, j)*»2+?.odoo*c< 5 , n* 

*C ( 4 , J ) *° ? ) 

J = l 

no TO N = 1 » 6 

S (1, J ) = E (1 • 1 )*XXINT (1 ,N )+E( l, 3)*( XYINT I l, NM-XYINT (N, 1) )+E( 3,3)* 

*yvintu.n) 

J=JM 

S(lfJ)*E(l,?)*XY!NT(l,N)+E<lf3)*XXINT< 1 , N ) +E (3 , 2 ) *YY INT ( l , N ) +E l 3, 3 

* ) *XY I NT ( N , l ) 

S(2,J)=E(2»ZJ*YYINT(l,N)+E(2,3)*<XYlNT<N,mXYINT(l,N))+E(3,3)* 
*XXINT(l, N) 

10 J = J + 1 
J = ^ 

HO 20 N=2,6 

S ( 2» J) = F (2,1 )*XYtNT (N, l ) + E( 2»3) *YY INTI 1»N >*F( 3. I >*XXINT( l,NH- 
*F( ^)*XYINT(l,N) 

S (3, J) = F( l, l)*XXINT(2.N)«-E( l,3)*< XY INT ( 2,N ) +XYI NT ( N, 2 ) )+E(3,3)* 

* YY INT ( 2, N ) 

J=J*l 

S(3t J ) =E { l ,2 )*XY!NT (2, N l ,3)* XX I N T ( 2*N ) +0 (2 : ,3 ) *YYINT( 2, N) + 

*F (3f 3 )*XYINT ( N , ? ) 

S<4,J)=F(2»?)*YYINT(2,N)4-E(2,3»*(XYTNT(N,2M-XYINT(2*N))*E(3,3)# 

* XX I NT (2»N) 

?0 J=J+! 

J = 5 

OC 30 N=3,6 

F (4, J) =F( 2, l)*XY INTIN, ? )>E( 2,3) *YYINT( 2,N>+E ( 3 , 1 ) *XX INT ( 2 , N ) + 

*F { 3 , * )*XYINT (? « N I 

M3,J)=FIl,l)*XXTNT<3,N)+E(l,3)’MXYINT(3,N>+XYINT(N,3)H- 

*F(3,3)*YYINT(3.N) 

J = J + 1 

C ( 5 » J ) = F { 1 * 2 ) * X Y I NT ( ^ » N ) + E ( 1,3)*XXINT(3,M)+F(2.2)*YYINT(3,N)«- 
* r (3, 3)*XYTNT(N,^) 

S(6, J) = E ( 2, ? )*YYINT(3,N) + E( 2,3)*(XYINT(N,3)*XYINT(3,N))+ 
*F(3,3)*XXINT(^,N) 
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j=j+] 

J = 7 

90 40 N=4,6 

S ( 6 , J) = R ( 2* 1 )*XYINT (N,3)+E(2,3)*YYINK3,N)+E(3,l)*XXINT(3,N) + 
( 3 » 3 )*XYINT(3,NI 

S (7, J) = E(lt n*XXINT(4,N) + E(l,3)*<XYlNT(4,Ni+XYINHN,4) ) + 

*E(3, 3l*YYINT(4 t N ) 

J=J+l 

S(7,J)=F(l,2>*XYINT(4,NM-E(l,3)*XXLM(4,N) + E(3,2)*YYINT(4,N) + 

* E ( 3 , 7)*XYINT(N,4> 

S(8,J)=F(2,?)*YYTNK4,N)*E(2,3)*(XYINT(N,4H-XYINT(4,N)H- 
*F( 3, 3 >*XXINT<4,NF 
40 J=J + l 
J =9 

nc 50 N= 5 » 6 

S ( 3 • J)=E(2.l )*XYINT(N,4)«-E(2,3)*YYINT(«,N)+E(3.l ) * XX INT ( 4 , N ) ♦ 

* E { 3 . 3 1*X Y I NT ( 4 , N ) 

S(<>.J)® f: (l,n , « I XXINT<5,N) + E(l,3)*<XYINT<5,NF+XYtNT<N,5)) + 
*E(3,3F*YVINK5,N) 

J=J + 1. 

St9,J)=F(l,2) *XY INT( 5,N )+E(l ,31 * X X I NT ( 5 , N } * F ( 3 , 2 > * YY INTI 5 ♦ N ) + 
*F(3»3)*XYINT ( N , 5 ) 

${10.J)«F(?,2)*YYINT(5,N)+E(2,3)*(XYINT(N,5) +XYINT (5 ,N > ) + 

*E (3, 3 ) *XXI NT ( 5 ,N ) 


50 J=J+l 

S ( 10 til ) =E< 2t l )*XYTNT (6. 5 ) + E( 2,3 )*YY INT (5 ,6 ) + F (3 ,1 »*XX INTI 5,6 ) + 
*E<3, 3)*XYINT( 5, A ) 

5 ( 11 , m=E( I, 1 ) * X X I NT ( 6, 6 ) + E( 1,3) *( XYINT( 6 ,6 U *2. ODOO+E ( 3 ,-3 ) * 
*YY INK 6, 6 ) 

SfU,l2»*E(l.2)*XYlNT(6,6) + E(l,3)*XXINK6,6)+E(2,3)*YYINT 
* ( 6 ♦ 6 ) + F ( 3,3 )*XYINT(f, 6) 

M 12, 1 21 = F (2, 2 ) *YYTNT { 6,6 H-E( 2,3 ) *< XYIM (6 ,6 )) *2.0D0Q+E< 3,* )* 
*XXINT (6, 6 I 
on 60 i = i,n 
K = I+l 


OC 60 J= K » 1. 2 
S(J, !)=$(! ,J) 


60 CDNTTNUF 

f, F.EE CeSS&Cf;C&f.E 6 f. 6 SF J GgEGG 6 r j 66 &C&C& 6 fl 666 SESS EEGEEeESE&SEEEE&EEE&eEE&EEEE 

■nr 70 1=1,12 


on 70 j= l , 12 

S( F, J) = S( I, J )*THICK 
70 CCNTINUE 

CESESEGBEE SB EE SEE S £ E E E E S E E E E E S E EE EEEE E EE £ EEE EEES EE S E E E E E EE EEEEEE BEE EE EEE 
PETURM 


FNO 
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SU^PHUTINP STRA in ( XG, YG,NQO,NELH ,NNP ,F , I 6 ,F f NRCW , STRSX, STPSY, STRS 
*XY,RTRNX ,STRNY, STRNXY J 

IMPLIfTT R F AL * R ( A-H, C~ 7. ) 

n IM P NS ION XGIMNP), YG(NNP) , NOD (NFL E ,6 ) , E ( 3 , T ) ,F ( NROW) »r(6,M 
niMFIMSION STRSX(NELF,6 ) , STRS Y <NF LF ,6 ) , ST RS XY ( NFIE, 6 ) ,STRNX « NFL E , h ) 
* , STRNY(MFL r , <: ) » STRNXY ( NELE » 6 ) 
niMFNSinN X(M *Y(6)tITl(6>. IT2(6> 

KNl=NnO( IP, l ) 

KN?=\ , n r '( r. c , .? } 

KNR=NCin ( m,*m 
KN4=NP0( if, A ) 

KN5=NOO( JF,S) 

KN6=NOO ( IF ,6 ) 
c $TABL IS h C MATR IX 
X (1) = XG( KM l ) 

X l ?. V=XG( KM’ ) 

V (1 ) = XG( KNT ) 

X (4) = XG ( KNM 
X (5I = XF( KMS 5 
X(6)=XG(KMM 

V { l) = YG( *Nl ) 

Y { ? ) = YG< KM’ ) 

YP )=YG( KMT) 

V ( a )~ yn ( kma. ) 

YM)=Y0( KMP ) 

Y(6)=Y T( KM *) ) 

rcMouT^ c matrix 
nr i i = i ,6 
r ( i , i ) = i .onoo 
cm ,’ ) = x i M 
c u , t ) - v ( m 
C t I , A ) = X f t)*Y( l ) 

ci i, si=x( n*x( t ) 

I C ( 1,6) =Y ( I I * Y ( r ) 

INVFPT C MATRIX AMO CVFR^RITE 
CALL 0MTNV(C,6,0,m, IT?) 

COMP'JTF STRAINS AT NO 06 S 
HC ’ F M=I ,6 
STRNXI TF,M)=n.OOCO 
STRNY UF,N) =O.OOCO 
STRMX Y ( I E ,N ) =0 , 0000 
nr ’S j= 1 , a 
N N =N n D ( I F, J ) 

J lj=? *NN~ 1 
JV=JLH-l 

STRMX ( 16 ,M ) = STRNXI TE, N ) -» ( C ( 2, J) +Y(N) *C< A, J) + 

*'’ o 0000*X (M I *C (Ft J ) J *F ( JU 1 
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STRNYUr,N)=STRNVUF,N) HC< 3 , J ) +X ( N ) +C ( 4 , J H 
«?.0 r »C0*Y(N) *C {4, J) )*F ( JV) 

^ STRNXYt TF,N) = STPNXYUF.*N) + (Ct 3, J) +X( N)*C ( 4« J )+2. 0D00*Y( N)* 

*f( 6 t J ) ) *F ( J- 1) + ( c ( ? , J H- Y ( N ) *C ( 4 T J ) +? . OD CO* X { N > *C ( 5 . J ) ) * F ( JV) 

CCMPUTF STRESSES 

nn ?. e . n=i»a 

STRSX HP ,N)=FI i t l ) *ST RNX ( I£»N) + E( 1 1 2)*STRNY( IF f N ) + E ( 1 , 3) * STRNXY ( 
* TF,N ) 

STR S Y ( I F » N ) =F ( i,;i)«STRNX( I F , N ) +E ( 2 . 2 )* STRNY ( IFtN)+E{2,3) 

**STRNXY ( l Ft N ) 

STRSXYnFtN)=F(l,3)*STRNX(IE»N)+E(2,3)*STRNY{ IF,NI4- 
*F(3t3)*STRNXY( IE.N) 

RET'JPN 

FNO 
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SLBRnUTINF SYRATU*B,NRCW t NCOL ,ISIG.NV,NC) 
r 

l Z TH I ^ SUBPD'J TINE SOLVES SYSTEMS OF SIMULTANEOUS LINEAR EQ- 

r. MATICNS FOP WHICH t>f MATRIX OF COEFFICIENTS (CALLEC A) 

C T S BOTH SYMMETRIC ANn BANDED. THE MATRIX A IS DIMENSION- 
r ED A (NROW ,M COL ) WHERE NP OW IS THE NUMRER GF EQUATIONS AND 

r. NCHL IS THF BANC WICTF, THE PRINCIPAL CIAGONAL IS STORED 
r in TEE F I D S T COLUMN OF A AND THF NCN ZERO BANDS ABOVE THE 
r PRINCIPAL DIAGONAL FORM THE BALANCE OF A. THE RIGHT HAND 
C SinE COLUMN VECTORS ARE STORED IN R WITH EACH P IGHT HAND 
C S IDF FORMING ONE COLUMN. HENCE E IS D T MFN S I ON c D BINROW, 

C NRHS ) WHFRF NR H S IS THF NUMBER OF RIGHT HAND SIDES FOR 
r WHICH THF SCLUTICN IS DESIRFD. KO IS AN EXECUTION INDIC- 
" ATPR. IF, UPON RETURN TO THE CALLING PROGRAM, KC IS ZERO, 
r T H c EXECUTION WAS SUCCESSFUL AND THE ANSWERS APF STOREC 
C IN R. IF KC IS NOT ZERO, THE EXECLTICN FAILED BECAUSE 
0 THF PRINCIPAL DIAGONAL CLEMENT IN THE KC’TH ROW WAS ZERO. 

C IN THIS EVENT \ NOTE IS WRITTEN CN THE OUTPUT PAGE. THIS 
r RPIJTINC DFSTROYS THE MATRICES A AND 8. 

IMPLICIT REALMS! A-H,C-Z ) 

DIMENSION A ( N V , N C ) , B ( N V » I) 

C MV, NO = ABSOLUTE DIMENSIONS CF A , B NROW = NO OF VARIABLES IN MATRIX 
0 NO PL = BANDWIDTH OF PRESENT SET OF EQUATIONS 

0 ISTG = ? SIGNALS ONLY &FDJCTIPN OF P.H.S. VECTOR AND BACK SUBSTITUTION 

101 FORMAT (//64H YOU GOOFED THERE IS A ZERO ON THF PRINCIPAL Q I AGO 

INAL IN TFE,l4 t 8H TH RCW.//) 

GO TO I IB , ID) , IS IG 
IB CONTINUE 

NHALF=NCCL-l 
NROWM=NR CW- l 
NRPW H=NR EW-NHALP 
) C ft 1 = 1 , NRO WM 
I F ( A ( 1,1)) I, IB, I 
1 P EC TP = 1 . OOOQ/A ( I ,1 ) 

IF(T-NROWH) 

? LIMIT = NH AL F 
GO TO A- 

B L I M 1 T=NP OW- I 
^ DC ft J =1 , L I M I T 
JROW= IfJ 

R AT [ P =- A ( I, J Fl ) *RFC ID 
i DO 3 K = J, LIMIT 

JCCL = KM -J 

B A ( J P OW , J CO L 5 = A ( JP OW , JCO L ) +R AT I G* A ( I , Kf l ) 

( 6 CONTINUE 

l.o CONTINUE 

c 
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HO 17 1=1 , NROWM 
no 17 J=l. .NHALF 
T J = I 4- J 

T F { IJ.OT .NRQW ) GO TO 17 

RUJ,U = 3 ( T J » l ) - R(|, i)*A( I*J + 1)/A( 1,1) 
17 CCNTINIJF 
r 

[F ( A (NROW T 1. ) ) 7, IA, 7 

7 RFCI D =l* CDOO/A(NROV»*l) 

3 R (NPOW, 1)=B( NROW, 1 )*PEC I P 

00 13 I=l,NRCW* 

1 R QW=NRO V> - 1 

«ECIP=1. OOOO/A { IROW,l ) 
t F ( [-NHALF) 0 , 10, 10 
n l I M I T = I 
GO TO 11 

10 l_IMIT=NHALF 

11 OC 1.2 J = 1 t L I w I T 
J RCW= IROW+ J 

J COL = J*l 

IP B ( IF CW ,1 )=6{ IRCW,l)-A( IROW, JC CL) *B (J ROW , l 1 
n P( TROW, l >=B( IRCW, 1)*PFCIP 
K 0 = 0 
P FT'IRN 
1A r=NPrw 
15 K C= I 

W° IT E ( 6* 101. )K0 

RETURN 

FND 


